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1. Introduction 

Parton Shower based Monte Carlo Event Generators (PSEGs) have developed into essential 
tools in High Energy Physics. Without them it is questionable if it at all would be possible 
to embark on large-scale experiments such as the LHC. Although they are based on leading 
logarithmic approximations and phenomenological hadronization models, they are typically 
able to describe hadronic final states in great detail and, especially at LEP, with great 
precision. However, there are problems. The description of final states which include 
more than three hard jets is not very good, and when it comes to collisions with incoming 
hadrons, the precision is generally lacking, especially for small-x processes. In this article 
we will address both these problems. 

The problem with describing several hard, well separated jets is inherent in the lead¬ 
ing log approximations, since they assume that there is strong ordering between parton 
emissions, and hence only give a good description of soft inter-jet and collinear intra-jet 
emissions. Typically it is possible to correctly describe one additional hard jet on top of 
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the hard sub-process used as starting point, by applying correction factors to the basic 
splitting functions. However, to go beyond one additional jet is more difficult. 

To describe events with several hard partons we can use so-called Matrix Element 
Generators (MEGs), where the parton distributions can be generated according to exact 
tree-level matrix elements. Unfortunately these matrix elements are divergent in the soft 
and collinear limits and a cutoff is needed to avoid these regions of phase space. However, 
to generate realistic events we need to hadronize the partons into jets of hadrons, and all 
reasonable hadronization models require that also soft and collinear parton emissions are 
modeled correctly. Hence the need for combining matrix element generators with parton 
showers. 

Gombining these two approaches is, however, not trivial. The matrix elements describe 
inclusive events, ie. events with at least n partons above some cutoff, while parton showers 
are exclusive and describes events with exactly n partons. Naively adding parton showers 
to events generated by a MEG will therefore give a very strong dependence on the cutoff 
used in the MEG, even if great care is taken to avoid double-counting by only adding 
parton showers below that cutoff. 

A solution for this problem was presented by Gatani et al. in [1]. The procedure, 
generally referred to as CKKW, relies on applying a jet clustering algorithm to partonic 
event from a MEG generating zero, one, two, etc. additional hard jets above some cutoff 
according to exact tree-level matrix elements. The repeated clustering of two jets into 
one is then used to construct an ordered set of scales corresponding to consecutive parton 
emissions. These scales are used to calculate Sudakov form factors corresponding to no¬ 
emission probabilities, which are used to reweight the MEG events to make them exclusive. 
A parton shower can then be added with a special veto to avoid double-counting of emissions 
above the cutoff. In this way it was shown that the dependence on the cutoff cancels to 
next-to-leading order accuracy. The dependence was, however, still quite visible, giving 
rise to annoying discontinuities in some observables. 

The basic CKKW prescription was improved [2] when implemented for the Dipole 
Cascade model [3,4] in the Ariadne program [5]. Rather than using a jet clustering algo¬ 
rithm to construct a set of scales, the Ariadne procedure involves constructing complete 
intermediate states corresponding to a series of emissions which the dipole cascade could 
have used to produce a given state obtained from a MEG. A special veto algorithm is 
used to calculate exactly the Sudakov form factors the cascade would have used to pro¬ 
duce the state, which are then used for reweighting. Together with a special treatment of 
the MEG-produced states with highest multiplicity, this procedure basically removes any 
visible discontinuities due to the cutoff. 

Both these procedures were originally developed for e'’“e“. Recently there has been 
some developments in applying them to hadronic collisions, in particular for the W-l-jets 
process, by Krauss et al. [6-8] and Mrenna and Richardson [9]. An alternative procedure 
has also been developed by Mangano [10], which is similar in spirit to GKKW, but which 
has a simpler interface between the MEG and PSEG. This development is very important 
for the LHG, where W-l-jets is an important background for almost any signal of new 
physics. 
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In this article we describe the extension of the CKKW procedure for the dipole cascade 
in Ariadne to handle hadronic collisions, again concentrating on the W+jets process. The 
goal is to obtain a procedure which gives as small cutoff dependence as was achieved 
for e'’'e“. However, we also expect to see differences w.r.t. the procedures of Mrenna, 
Richardson and Krauss, since the dipole cascade model for collisions with incoming hadrons 
[11,12] is different. The standard initial-state parton shower approaches, such as those 
implemented in PythiA [13,14], HERWIG [15] and Sherpa/APACIC++ [16,17], as well as 
the Sudakov form factors used in CKKW, corresponds to a DGLAP-resummation [18-21] 
of leading logarithms of the hard scale. In the dipole cascade, however, also some terms 
corresponding to logarithms of 1/x are resummed. Although not formally equivalent to 
neither BFKL [22-24] or CCFM [25-28] evolution, it has proven to be able to describe 
most features small-x final states at HERA, where all DGLAP-based parton showers fail. 
Now, W-|-jets is conventionally not considered to be a small-x process because of the hard 
scale, mW) being large, but at the LHC the collision energies are so large there may be 
substantial effects of terms proportional to a” \og{^fS/, where rriT^/^/S ~ x ~ 0.005. 

The layout of this article is as follows. We will first recap the main points of the 
CKKW procedure and how it is implemented in Ariadne for e'’'e“ in section followed 
by a description of how hadronic collisions are treated in Ariadne in section Then, 
in section we will describe our CKKW implementation for W-l-jets, and hadronic colli¬ 
sions in general, starting with the construction of intermediate states in O, followed by 
the reweighting procedure in O. The results of our investigation of its performance are 
presented in Finally in section ^ we present our conclusions. 


2. CKKW and the Dipole Cascade Model 

When generating events with a PSEG, the procedure is to start from a primary hard sub¬ 
process, typically a 2 —> 2 process such as e'''e“ —> qq or qq —> W —> ue, and then to let the 
incoming and outgoing quarks and gluons evolve a parton cascade in an iterative 1 —> 2 
branching procedure. The emissions are ordered according to some evolution scale p, where 
the maximum scale, po is typically given by the hardest scale in the primary sub-process, 
and the minimum is some cutoff scale of the order of one GeV, typically tuned to match a 
particular hadronization model. 

We can write the exclusive cross sections for in this way generating 0,1,2,... additional 
partons above the cutoff, pc, as 

o'+o = o-Q ^SoiPO, Pc) 

d(7+i = ao as{pi)cif Asoipo, pi)/^Si{pi, Pc)dpidQi 

da+2 = aoas{pi)a^{p2)c\2 ^SoiPo, Pi)^sApii P2)^S2{P2, Pc)dpidVtidp2dVL2 

n 

da+ri = aoc^^As„{pn,Pc)Y\as{pi)Asi_i{pi-i,Pi)dpidQi 

i=l 

; ( 2 . 1 ) 
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where the ordering \s pQ > pi > ... > pn > Pc and fij symbolizes the phase space variables 
defining the ith emission in addition to pi (typically some momentum fraction, Zj, and 
some azimuth angle 4>i). As^{pi, Pi+i) is here the so-called Sudakov form factors giving the 
probability that no emissions occurred from the state with i additional partons between the 
scales Pi and Pi+i- The coefficients are basically products of splitting functions which 
depends on pi and Oj, and we assume an implicit sum over all possible flavour combinations. 
As we will see below, the may also include ratios of parton density functions (PDFs) 
in the case of incoming hadrons. 

The Sudakov form factors are formally resummations of virtual diagrams to all orders 
and can, in principle, be calculated analytically. They would then only depend on the 
limiting scales, as is done in the standard CKKW procedure. However, when explicitly 
interpreted as a no-emission probability in a PSEG, as is done eg. in Ariadne, it basically 
depends on all momenta in the partonic state. The typical form of the Sudakov is 

^s{Pi,Pi+i) = ew (^- J ^as{p) J dzP{z)^ , ( 2 . 2 ) 

which, of course, can be expanded in a series in Og, and we can rewrite the exclusive cross 
sections as 

<7-1-0 = <7o (1 + Cqi as + Cq 2 ag -h ...) 

dc-f-i = iToQ;scrf(l + cf;fas + T • • ■)dpidfli 

da+2 = o-Q al(^i{l + C23®as + Q^s + • • ■)dpidVtidp2dVt2 

n 

da-i-n <7o Og C^^(l -|- “k Cri,n+2*^s T . . .) dpid^i 

i=l 

; (2.3) 

to emphasize the resummation aspect. We note that even though all the coefficients c?® 
are divergent in the soft and collinear limit when pc ^ 0, the resummation to all orders in 
the Sudakovs gives a finite result for each of the cross sections. Also, when integrated over 
the allowed phase space the cross section of the primary sub process ao is retained, 

oo 

^fj+i = fTo. (2.4) 

0 

In contrast a MEG will generate inclusive partonic states with the cross sections for 
generating at least 0, 1, 2, ... additional jets given by 
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O'+O — O-Q 

da+i = (To asCn^(irii 
(i(T+2 = aoagC^dClidfl2 

n 

da+n = ctq U dfti 

i=l 

■■ (2.5) 

where the symbolizes all phase space variables defining the fth parton, and the coeffi¬ 
cients are calculated using the exact tree-level matrix elements (including PDFs in the 
case of incoming hadrons). Clearly we can not simply add these cross sections, especially 
since each of the coefficients are divergent if the soft and collinear limits are not cut off 
properly. 

The advantage of using a MEG is that the exact tree-level matrix elements are used, 
which means that also states with several hard partons are described correctly. This is not 
the case for the PSEG, where the coefficients are given by products of splitting functions, 
which is only a good approximation in the limit of strongly ordered emissions. On the 
other hand a MEG will not correctly treat soft and collinear partons, where the coefficients 
are large and need to be resummed to all orders, as in a PSEG, to give reliable results. 
Glearly it would be highly desirable to combine the two approaches. 

It should be noted that for the first emission in a PSEG is typically quite easy to 
modify the splitting functions to correctly reproduce the exact matrix element, effectively 
replacing cff with and in most PSEGs this is the default behavior for most primary 
sub-processes [3,5,12,29-37]. 


2.1 The original CKKW procedure 


Gomparing eqs. and the solution should be obvious. Use a MEG to generate up to 
N additional partons above some cutoff, r/cut, but reweight the generated states with the 
Sudakov form factors, and then add a parton shower with the requirement that no partons 
above Ucut are emitted. This is the essence of the CKKW procedure. To calculate the 
Sudakov form factors we need an ordered set of emission scales, which is not provided by 
the MEG, since there all possible diagrams are added coherently and emission scales are 
not well defined. In the original CKKW procedure, the A:_l clustering algorithm [38,39] was 
used to define an ordered set of scales which were used to analytically calculate the Sudakov 
form factors. A reweighting was also done to have the constructed scales as argument to as- 
The resolution variable of the A:_L-algorithm was also used for the cutoff in the MEG, which 
is not the same as the evolution variable in the PSEG. To ensure a full coverage of the phase 
space the parton shower was therefore added with the maximum scale as starting point, 
but vetoing all emissions corresponding to the A:^-algorithm resolution variable above j/cut 
to avoid double-counting. 
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It was shown that this procedure removes the dependence on the MEG cutoff, ycut, to 
next-to-leading logarithmic accuracy. However, there was still a clearly visible discontinuity 
in some generated distributions. In [2] the procedure was improved in several ways when 
it was implemented for the AriADNE program. 


2.2 The Dipole Cascade Model 


Ariadne implements the dipole cascade model which is quite different from conventional 
parton cascades. Rather than iterating 1^2 parton splittings, gluons are emitted from 
colour-dipoles between colour-connected partons resulting in 2 ^ 3 parton splittings. This 
model has several advantages. Since gluons are emitted coherently by colour-connected 
partons, there is no need for explicit angular ordering. In addition, the evolution variable 
is defined as a Lorentz-invariant transverse momentum which also is a suitable scale to be 
used in as- The evolution variable is defined as (for massless partons) 


pI = 


■§12523 

5123 


( 2 . 6 ) 


where parton 2 is the emitted one and Sij and Sijk are the squared invariant masses of the 
two- and three-parton combinations. 

The probability for a given emission is given in terms of the dipole splitting functions, 
which depend on and a Lorentz invariant rapidity defined as 


1 , 512 

y = o In-. 

2 523 


(2.7) 


The probability of a gluon emission from a dipole between two partons i,j is then given by 


dP{p\,y) = as{p\)Dij{p\,y) exp 



dy'as{pl)Dij{pl,y' 


'2 


dp\ 

-^dy, 

PI 


( 2 . 8 ) 


where exp(...) is the Sudakov form factor. The dipole splitting functions, Dij, depends on 
which partons are involved according to 


-^qq(P±)y) ~ 


Xf -|- Xg 


a 


qg(7'±)y) ~ 


Pggip±iy) — 


Stt (1 - Xl)(l - Xs) 
3 Xf -|- Xg 

dvr (1 - xi)(l - Xs) 
3 Xf -|- Xg 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


dvr (1 - xi)(l - Xs) 

where x* are the resulting energy fractions of the emitting partons in the original dipole 
rest system, Xj = 2Ei /related to p^ and y according to 


y = pi = si 23 {l - xi){l - Xs). (2.12) 

It can be shown that the dipole splitting functions are equivalent to the standard Altarelli- 
Paris! splitting functions in the relevant soft and collinear limits [3]. We also note that the 
Dqq exactly corresponds to the leading order e"*“e“ ^ qgq tree-level matrix element. 
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Another feature of the dipole cascade, which will turn out to be important for the 
CKKW implementation, is that all partons are always on shell throughout the cascade. 
This is possible since the recoil from the emitted parton can be absorbed by the two emitting 
ones. In contrast, a conventional parton cascade does not have on-shell intermediate states, 
and the full kinematics of an event is not constructed until all the scales in the complete 
shower have been generated. While the energy loss of the emitting partons are defined in 
the splitting functions in eqs. (2.£)-(2.11), the transverse recoil is chosen according to some 
principles detailed in [5]. 

It can be noted that the “inverse” of the dipole cascade is a well-behaved jet clustering 
algorithm. In fact such an algorithm has been constructed, the DiCLUS algorithm [40], 
based on successive clusterings of three jets into two, using the p±m eq. (|2.6|) as resolution 
scale, which has been shown to have many attractive features [41]. 

There are, however, also some disadvantages with the dipole cascade model, the main 
one being that it only deals with gluon emissions, and the splitting of gluons into qq pairs 
must be added by hand, both for final-state [42] and initial-state [34] splitting. 

The initial-state splitting will be described in section |3.2| below. Final-state splittings 
are simply added as a possibility for a dipole connected to a gluon to split this gluon into a 
qq-pair in addition to emit a gluon. Here, the standard Altarelli-Parisi splitting function 
is used, divided between the two dipoles connected to the gluon. This will result in a new 
dipole splitting function, 


= 

*g 4vr 


^ (1 - X 2 )^ + (1 - 


1 — Xl 


(2.13) 


In the original formulation, the splitting was divided equal between the two dipoles con¬ 
nected to the gluon, ie. ^ = 0.5. However, in the current Ariadne implementation, a larger 
fraction is given to the smaller of the dipoles ig and gj, hence for the ig dipole we have 

^ ~ ^sj /+ '^*g)- 

Although these splitting do not come in naturally in the dipole picture, they can be 
incorporated in a consistent way and the resulting implementation in Ariadne is probably 
the best model for describing both e'’'e“ final states at LEP [43], and DIS final states at 
HERA [44]. 


2.3 Ariadne and CKKW 

The Ariadne implementation of CKKW for e'’'e“ is described in detail in [2]. Consider¬ 
ing the nature of the dipole cascade it may seem reasonable to use the dipole clustering 
algorithm in DiCLUS to construct scales, rather than using the A:_L-algorithm. However, in 
the Ariadne implementation a further step is taken. For each partonic state generated 
by a MEG, all possible dipole cascade histories are constructed, basically answering the 
question how would Ariadne have generated this state? A specific history is then picked 
by weighting possible histories with the product of the corresponding dipole splitting func¬ 
tions. The implementation depends on the MEG generating specific colour connections 
among the partons. Although this information is not physical, it is usually provided by 
ME generator programs. (See discussion in ref. [45]). In principle one could choose between 
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all possible colour connections in the same way as different histories are considered, but as 
most MEGs supply colour information this is not necessary. 

From a MEG generated state with n additional partons, we can now construct, not 
only an ordered set of emission scales, corresponding set 

of intermediate states, 5^-1, Sn-2-, ■ ■ ■ Sq. As in the standard GKKW, the reweighting with 
the correct scales in is done using just the constructed scales. However, the Sudakov 
form factors are calculated using the fact that exactly corresponds the 

probability that no emission occurred from state Si between the scales and 
Hence, letting Ariadne make a trial emission, starting from the state Si with p'^- as the 
maximum scale and throwing away the MEG event if the trial emission was above 
will exactly correspond to reweighting with the same Sudakov form factor, A^. 
which Ariadne would have used when generated the event. 

A special treatment is given the trial emission from the MEG generated state, Sn- 
Rather than using the cutoff, from the dipole cascade, the event is thrown if the 
emission is above the cutoff used in the MEG, while if the emission is below, the emission 
is kept and the cascade is continued down to p^^ to produce a full ME+PS event. In 
addition, if n = iV, the highest number of additional partons generated from the MEG, the 
emission from S^ is always kept. This was not done in the original CKKW prescription, 
but is clearly needed, since otherwise we would never get events with N + 1 additional 
partons above the cutoff. In later developments of CKKW such a treatment has been 
added [7,9,46]. 

The whole procedure looks as follows: 

1. First the number of partons, n < N, to be generated is chosen according to the 
integrated tree-level matrix elements in the MEG, using a cutoff ycut in some jet 
resolution scale and a fixed as(Pj_c)- 

2. Then MEG is told to generate the momenta of the state with n additional partons 
according to the tree-level matrix element. Since we do not want any events below 
the cutoff in the dipole cascade, the invariant of the partons is checked, and if 
anyone is below the state is rejected and the procedure is restarted at step 1. 

3. Now, all the intermediate states 5^-1, ..., 5o and scales p\^, ..., constructed 

corresponding to a possible dipole shower history of the generated Sn state. 

4. The generated event is rejected and we restart at step 1 with a probability given by 
Hi ®s(PAj)/os(p±c)”' to get the correct scales in ag 

5. We now make a trial emission with the dipole cascade from the state Sq starting 
from the maximum scale, typically limited by the squared center-of-mass energy. If 
this emission is at a scale above the event is rejected and we restart from step 
1. If not, a trial emission is performed from the state Si with a maximum scale of 
p^i- If this emission is at a scale above p‘j _2 the event is rejected and we restart from 
step 1. This procedure is repeated for all states down to 5^-1. If no rejection has 
been made, a trial emission is made from the ME-generated state with n additional 
partons starting from the scale There are now two cases: 



(a) If n = the trial emission is always kept and the dipole cascade is allowed to 
continue down to the cutoff and the event is accepted. 

(b) If n < and all parton pairs pass the cut, ycut used in the MEG, the event is 
rejected and we restart from step 1. If any of the partons fail the cut, the trial 
emission is accepted and the dipole cascade is allowed to continue down to the 
cutoff and the event is accepted. 

Taken together, this so-called Sudakov veto algorithm, giving “exact Ariadne Su- 
dakovs”, and the special treatment of the states with highest multiplicity, results in a 
dramatic improvement of the smoothness of observables sensitive to the cutoff region for 
e'''e“ events [2]. 


2.4 CKKW in hadronic collisions 

To extend this algorithm to be used for hadronic collisions is conceptually straight forward. 
In [7] and [9] the standard CKKW procedure was extended to hadronic collisions in general 
and for the W-l-jets process in particular, and implemented for the APACIC++, HERWIG and 
Pythia PSEGs. The principle is the same as for e'*'e“. Jet construction is done with the 
A;_L-algorithm modified for hadronic collisions according to [47]. The resulting ordered set 
of scales is used in the analytic Sudakov form factors (in [9] a Sudakov veto algorithm 
similar to the one in Ariadne was used for the PYTHIA implementation), and the events 
from the MEG was reweighted with these and the properly scaled ag- The MEG states 
with highest multiplicity was treated in the same way as in the Ariadne implementation 
above. 

There are some issues which need to be treated with special care. In a PSEG, the 
initial states emissions are generated in a backward evolution procedure which besides the 
standard partonic splitting functions also involve ratios of PDEs. The leading order cross 
section is given by 

dao = dahh^w = (2.14) 


where the scale fj? = = Xj^X-S. Making one step in the backward evolution with a 

g ^ qq splitting will be performed with a probability 


dP{Q^,z) = asPg^diz) 


X+fq{x+,Q^ 


:Aa(QLx,Q') 




dz, 


(2.15) 


where the Sudakov form factor can be formulated both with (Pythia) and without 
(herwig) ratios of PDEs. Both choices are formally equivalent in the leading-log ap¬ 
proximation, but only the former choice corresponds exactly to a no-emission probability. 
The maximum scale is typically given by and the Sudakov form factor corresponds to 
a leading-log DGLAP resummation^. However, clearly there is nothing in the real world 

^Also next-to-leading logarithmic Sudakov form factors may be used [1], although these may become 
larger than unity, disabling the interpretation as no-emission probabilities. 
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preventing emissions with a scale above and such parton states will be generated by 
the MEG. In [7] states with one or more partons emitted at a scale above are treated 
as coming from a separate class of primary sub-processes, and the Sudakov form factors 
are added only for additional emissions with a maximum scale given by the smallest of the 
constructed scales above 

Although the resulting procedures in [7] and [9] are shown to be fairly cutoff insensitive, 
there is still some dependence, and it may be worth while to extend also the Ariadne 
implementation of CKKW to hadronic collisions to see if a better result can be achieved. 
As for the standard CKKW this extension is in principle straight forward. However, there 
are some tricky issues, mainly to do with the treatment of parton densities, and to describe 
how we deal with these, we first have to describe how hadronic collisions are implemented 
in Ariadne. 


3. The Dipole Cascade Model for Incoming Hadrons 

In contrast to conventional parton shower models, the dipole cascade model does not 
separate between initial- and final-state gluon radiation. Instead, gluons are always emitted 
from final-state dipoles as in the e^e“ case. The cleanest situation is in the DIS electro¬ 
production case, where the leading order process is eq —> eq, ie. a quark is being kicked 
out of a hadron. A gluon may then be emitted from the colour-dipole between the struck 
quark and the hadron remnant, using the same dipole splitting function as in the e'’'e“ 
case in eq. ( [2.91 ). There is one major difference though. In the e'’'e“ case, both the quark 
and anti-quark can be considered point-like, but in the DIS case only the struck quark is 
point-like (at least up to the resolution scale, of the exchanged virtual photon) while 
the hadron remnant is an extended object with a size of roughly one fermi. Just as for the 
electromagnetic case, radiation of wavelengths much smaller than the size of the antenna 
is suppressed. 


3.1 Gluon emission from an extended source 

In [11] it was argued that only a fraction of the hadron remnant is effectively taking part 
in the emission. For a gluon emission at the scale p\ this fraction is given by 


a{p±) = 



(3.1) 


where p, parameterizes the inverse size of the remnant and a reflects the dimensionality of 
the emitter. If only that fraction of the remnant momentum is allowed to take part in the 
emissions, this corresponds to a sharp cut in the allowed phase space for gluon emission, 
limiting the transverse momentum mainly in the remnant direction according to 


P± < 


Wajp^) 

e~^y -|- a{p±)e~y 


(3.2) 
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Here we note another major difference as compared to conventional initial-state parton 
showers. From eq. (3.2) it is easy to see that the maximum scale is given by 


P±ma,x — 




(3.3) 


which may be much larger than which is used as the maximum scale in conventional 
PSEGs, especially for small x values since jx. Since the maximum scale is also 

used in the Sudakov form factors, these do not only correspond to a standard DGLAP 
resummation of leading logarithms of but in the dipole model they also resum, at least 
partially, logarithms of 1/x. Note however, that there is no formal equivalence to BFKL 
or GGFM evolution. There is another similarity though. Even though the emissions in the 
dipole model are ordered in p±, they are not ordered in rapidity, and conversely, following 
the emissions in rapidity from the struck quark the transverse momenta of the gluons will 
be unordered as in BEKL and CCEM. 

In the implementation in Ariadne, the sharp cutoff in p±is replaced by a smooth 
function, 0(p^,?/), with some power suppression for emissions violating eq. ( |3.2|) . Goncen- 
trating on the soft and collinear limits, where the splitting function is simply oc dlnp‘j_dy, 
we can write the probability of emitting a gluon as 

dP{p‘i,y) = ^Q{p'i,y)As{W‘^,p‘i)^dy (3.4) 

37r p. 


Comparing this to the corresponding initial-state q —> q splitting in a conventional parton 
shower, (cf. eq. (|2.15|)) where we have 


dPiQ^,z) = 


4as 1 ^/q(^,Q"), 


37r 1 -2;x+/q(x+,Q2) 


^s(f5max>Q 




dz, 


and noting that in this limit 


1 


dQ^ dp\ 

-dz = — 7^dy, 


p\ 


z{l-z) Q2 

we see that the suppression function, 0, corresponds to the ratio of PDEs 


0(p±,y) ^ ^ 


x+fci[x+,Q‘^) 


(3.5) 


(3.6) 


(3.7) 


The fact that only a part of the remnant takes part in a gluon radiation also means 
that only a fraction of it will obtain a transverse recoil in an emission. This is handled by 
the addition of so-called recoil gluons and is described in some detail in [5]. These recoil 
gluons will not be relevant for this report, however they will play a role when implementing 
CKKW in Ariadne for DIS and we will come back to them in more detail in a future 
publication [48]. 

For W production in hadronic collisions, the primary sub-process is qq —> W and the 
initial dipole from which gluons are radiated is between the two remnants. The model is 
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now the same as in DIS. However, since both remnants are extended, the cutoff in eq. ( |3.2| ) 
becomes 


P± < 


Wai{p±)a2{p±) 
a2ip±)e+y + ai{p±)e-y ’ 


(3.8) 


and the maximum scale is 


P±ma.x — 




1 

2+ai+a2 


(3.9) 


Again the sharp cutoff is replaced by a power suppression of transverse momenta above the 
limit in eq. (p.8|). Rather than introducing recoil gluons to absorb the transverse recoil, we 
note that any emission from the primary dipole corresponds to initial-state radiation, for 
which it is natural that the recoil is taken by the W (otherwise it would not be possible 
to produce a W with non-zero transverse momentum). In [12] the choice was to always 
transfer the transverse recoil to the W in the first emission, while in subsequent emissions, 
the recoil is only transferred if the emitted gluon is close to the W in phase space. 


3.2 Sea-quark emissions from remnant dipoles 


Besides gluon radiation, there is also a possibility that one of the quarks fusing into the W 
is a sea-quark, in which case it could have come from a perturbative splitting of a gluon 
as in eq. ( 2.15| ). As for the final-state gluon splitting into qq, this process does not come 
in naturally in the dipole model. Instead it is added by hand as an explicit initial-state 
splitting. This procedure is detailed in [34], and is based on different treatments of the 
remnants depending on whether a valence- or a sea-quark entered into the primary sub¬ 
process. A sea-quark is picked with the probability xfsq{x, Q'^)/{xfsq{x, Q^) +xfyq{x, Q^)) 
and in this case the complex remnant containing the anti-sea-quark and the valence quarks 
is split into a colour-singlet hadron containing the anti-sea-quark and a simple remnant. 
The sharing of longitudinal momentum is inspired by the string fragmentation function as 
explained in detail in [11]. A dipole connected to such a remnant is now allowed to emit 
the anti-sea-quark in a way similar to a standard initial-state parton shower, except that 
the ordering is in transverse momentum, changing eq. (2.15) to 


dP{p] 


— UoRp 


g^q 


£+/•('£+ „2 

z Jsy z A 2', 

(^) „ j; _2 \^^^P-i-raax^P±) 


X+fq{x+,Pi) 


dpi 

pI_ 


■dz. 


(3.10) 


As for the case of final-state gluon splitting, there is now in the W-production case sev¬ 
eral competing processes which can occur in the primary dipole, and after generating one 
emission of each, the one which gave the largest p± is chosen. If the emission of an anti¬ 
sea-quark is chosen, it will form a new dipole with the remnant of the hadron to which it 
previously belonged. The transverse recoil is taken by hard subsystem, just as in a stan¬ 
dard parton shower, where the hard subsystem in our case is the W and any other parton 
which has been previously emitted. 

It may seem counterintuitive that the essentially non-perturbative splitting of the rem¬ 
nant is allowed to affect perturbative emissions. However, the splitting will mainly influence 
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the region very close to the remnant, which is typically out of reach for current experiments. 
Nevertheless, we will investigate alternative treatments in a future publication. 

Emitting an anti-sea-quark means that we have now extracted a gluon from the in¬ 
coming hadron. In a standard parton shower scenario it would then be possible to evolve 
this gluon backwards either with a g —> g splitting or a q —> g splitting. In Ariadne, the 
former is modeled by gluon emissions from either of the two dipoles connected to the two 
remnants. In this case we use the same suppression function as in eq. (3.4), but comparing 
to eq. (3.10|) we find that it now corresponds to the ratio of gluon densities. 


Q{p‘i,y) 


x+fgix+,Q‘^)' 


(3.11) 


where the extra factor z in eq. (^^) is absent since this is now included in the gluon splitting 
function. The initial-state q —> g splitting is not included in the Ariadne program, but 
it could in principle be added in the same way as the sea-quark emission. Again we will 
investigate this in a future publication. 

Clearly, the dipole model for incoming hadrons has some conceptual problems, espe¬ 
cially when it comes to initial-state g —> q splittings. However, it also has some advantages. 
The hrst emission is quite easily modified to correctly reproduce the leading order matrix 
element, both for DIS and W-production. Also, a larger part of phase space is available for 
gluon emissions as compared to DGLAP based initial-state parton showers. This enables 
Ariadne to reproduce small-x observables in DIS, such as the forward jet rates, where 
no conventional shower succeeds. Below we shall also see that Ariadne gives a somewhat 
harder peak in the W p_L-spectrum than conventional parton showers, which we know peaks 
below the data measured at the Tevatron. 

Now that we have explained how Ariadne handles W-production in hadronic collisions, 
we can proceed with describing how to combine it with a fixed-order tree-level MEG. As 
for the e'’'e“ case it will involve constructing all possible cascade histories of a produced 
MEG state, the reweighting with Sudakov form factors using a Sudakov-veto algorithm, 
and finally the reweighting with as as well as with ratios of PDEs and the suppression 
function, 0. 


4. Ariadne and CKKW for W production 

The states delivered by a MEG contains information about the momenta, colour connec¬ 
tions and types of the incoming and outgoing particles in the generated sub-process. The 
states are generated according to exact tree-level matrix elements, using a fixed Og eval¬ 
uated at some scale Qq and weighted by the relevant parton densities typically evaluated 
at the same scale. Qq is usually taken to be the cutoff scale used to regularize soft and 
collinear divergencies. 

4.1 Constructing the Emissions 

To construct a dipole cascade history of this state it is first necessary to introduce the 
remnants so that all outgoing coloured particles from the sub-process are connected with 
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Figure 1: Different ways dipoles are connected (dashed arcs) depending on which kind of parton 
is extracted from a baryon. Filled ovals corresponds to colour-singlet hadrons, while open ovals 
represents coloured remnants. 


dipoles. This is done in the same way as for standard Ariadne. In figure |T| the different pos¬ 
sible connections of dipoles to the remnants are described schematically for W-production 
in pp collisions. We see that if a gluon has been extracted from a baryon (lower part 
of figure |I]b), there are two remnants one containing a quark connected to parton in the 
end of anti-colour line of the gluon, and one di-quark connected to the end of the colour 
line. Furthermore, if a sea-quark is extracted from a proton (upper part of figure ||a), the 
parton in the end of its colour line will be connected with a di-quark remnant, while the 
anti-sea-quark will form a hadron together with the remaining valence flavour, as described 
in section |3.2| above. Similarly, if an anti-sea-quark is extracted (upper part of figure mb), 
the parton in the end of its colour line will be connected with a single quark remnant, while 
the sea-quark will form a hadron together with the remaining valence flavours. Finally if 
a valence quark is extracted (lower part of figure ||a), the remnant is a di-quark which is 
connected with the parton on the end of the colour line. 

The construction will now proceed iteratively, each step corresponding to the inverse of 
an emission in the dipole cascade. All possible constructions will be made, and afterwards 
one of them will be picked. In each of the construction steps we must determine 

• The scale of the corresponding emission. 

• The value of splitting function to be used to give different weights to different possible 
construction paths. 

• The ratio of PDFs or the value of the suppression function which would have been 
used in the corresponding emission. This will be used to reweight the events. 

• The way the momentum of the emitted parton is distributed among the emitters. 
There will usually be three partons constructed into two, which means that the total 
energy and momentum is always conserved with all partons staying on-shell. However, 
the orientation of the final two partons in the rest system of the construction needs 
to be specified. 

Some of the possible construction steps correspond to emissions from dipoles between 


- 14 - 


















partons from the hard sub-process, and these are the same as in e'’'e“. Then there is 
a group of construction steps which involve hadron remnants, which are particular to 
hadronic collisions in general and to W-production in particular. In appendix^ we present 
a complete list of all possible construction steps 

After the construction procedure we are normally left with several possible cascade 
histories. Most of these will end up in a zeroth order state containing only remnants and 
one W with no transverse momentum. There will be some diagrams generated by the MEG 
which never could have come from an Ariadne cascade. One example is the initial-state 
q —> g splitting discussed in the end of section in which case the constructed state is 
accepted anyway. However, there are also diagrams, such as the one in figure which could 
not be produced even by a conventional parton shower. In this case the construction is 
stopped before reaching the zeroth order state, and this state is then treated as a separate 
leading order process and the reweighting is only applied to additional partons (similarly 
to the treatment in [7] mentioned above). 

The resulting alternative cascade 
histories may or may not have an or¬ 
dered set of constructed scales. When 
choosing a history according to their 
weights given by the products of the 
splitting functions, we first only con¬ 
sider true Ariadne histories with or¬ 
dered scales. Only if no such histories 
were found, the other histories are con¬ 
sidered. Histories corresponding to fig¬ 
ure ^ will only be considered if no full 
constructions are found. 



Figure 2: An example of a W-strahlung dia¬ 
gram. Such diagrams are not modeled by standard 
Ariadne. 


4.2 Reweighting the Events 


For a given MEG state, Sn, with n additional jets, we have now constructed a dipole cas¬ 
cade history with complete intermediate states, Sn, ■ ■ ■, Sq, and the corresponding emission 
scales, p±n, ■ ■ ■ ,P±i, and we can proceed with the reweighting. 

First we note that the MEG has used PDFs typically evaluated at the cutoff scale, Qq 
with and x_ given by the light-cone momentum fractions of the partons, i and j, entering 
the hard sub-process. This should be compared with the starting point for a normal parton 
cascade generation, where we just have a qq' —> W sub-process, and the PDFs are evaluated 
at the scale and and x'_ given by the corresponding momentum fraction for the q 
and q'. Our strategy is to follow the Ariadne cascade as closely as possible, just replacing 
the product of dipole splitting function with the exact tree-level matrix element, so to get 
the same starting point, we take the q and q' of the state Sq and their x'^ and x'_, and 
reweight the event with 


• x'_/q/ {x '_, m^) 
x+fi{x+, Ql) ■ x-fj{x-, Ql) 
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If a construction instead ended in a state such as the one in figure the corresponding 
incoming partons and their momentum fractions are used instead with the scale given by 
the squared invariant mass of the hard sub-process. 

Then we reweight with all the PDF ratios, determined in the construction, 

n 

= (4.2) 

i=l 


This comes about since the exact tree-level ME used corresponds to the product of splitting 
functions, while in a parton cascade we also have ratios of PDFs as in eqs. ( |3.5D and ( 3.1C| ). 
Depending on the emission, these ratios can be either 1 for a final-state emission, the ratio 
of PDFs for the case of initial-state q —> g and g —> q splittings, the suppression function 
0 for an initial-state g —> g splitting and Q/z for an initial-state q —> q splitting (cf. eqs. 
dO) and &)■ We note that for a conventional parton cascade, where the 0 functions 
would be replaced by ratios of PDFs, the loq and cui weights would basically cancel each 
other, which is why these did not show up in the procedures in [7] and [9]. 

We then reweight with the correct scales in cts according to 


W2 — 


nr=i«s(pL) 


(4.3) 


«s(Qg)- 

Again, for the situation in figure |2|, the first two scales are taken to be rrijj. 

Finally we need to reweight with the Sudakov form factors in Ariadne. This is done 
with the same Sudakov-veto algorithm as was presented in section 2.3. There are, however 
a few details which should be mentioned. 

The starting scale for the trial emission from the leading order state, Sq, is given by 
Plmax ~ IF^/4, where W is the total invariant mass of the hadronic collision, ie. the same 
as for the standard Ariadne treatment of W production. For the situation in figure rnjj 
is used instead. 

If the constructed cascade history contains unordered scales, such that 
the two corresponding emissions will be treated as a combined emission with as the 

scale. The Sudakovs will be generated with a trial emission from the state Sj_i with a 
minimum scale of and a trial emission from the state S'j+i with a maximum scale of 

Pii+i’ '''^hile there is no Sudakov generated from the state Si. 

Finally, in the trial emission from the state Sn, for n < N, when checking if the 
resulting partons are above the jet cutoff used in the MEG, possible recoil gluons are not 
considered. Such recoil gluons may appear in Ariadne, but they are typically rather soft, 
and including them would very often result in the emission being below the cutoff, even if 
the emitted gluon is not. 


4.3 The Full Algorithm 

We can now summarize the full algorithm. The way it is used below results in weighted 
events. This is because of the complicated reweightings which takes place. However, all 
weights are positive, and by carefully choosing the PDFs and Og used in the MEG, it should 
be possible to have a vetoing procedure so that all events end up with unit weight. 
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1. First the number of partons, n < N, to be generated is chosen according to the inte¬ 
grated tree-level matrix elements in the MEG, using a cutoff Qq in the jet resolution 
scale given by the longitudinally invariant /c^-algorithm. A fixed ag is used and the 
PDFs are typically sampled at the Qq scale. 


Then the MEG is told to generate the momenta of the state with n additional partons 
according to the tree-level matrix element. Since we do not want any events below 
the cutoff in the dipole cascade, the invariant of the partons is checked, and if 
anyone is below the state is rejected and the procedure is restarted at step 1. 


3. Now, all the intermediate states Sn-i, ■ ■ ■, Sq and scales 


, p^^ are constructed 


according to the procedure in section 4.1, resulting in a possible dipole shower history 
of the generated state. 


4. 

5. 


The event is reweighted by the weight factors given in eqs. (4.1)-(4.3). 


We now make a trial emission with the dipole cascade from the state So, starting 
from the maximum scale = lT^/4. If this emission is at a scale above p^i, 

the event is rejected and we restart from step 1. If not, a trial emission is performed 
from the state with a maximum scale of Pj^^- If this emission is at a scale above 
p ‘]_2 the event is rejected and we restart from step 1. This procedure is repeated for 
all states down to Sn-i- If no rejection has been made, a trial emission is made from 
the ME-generated state with n additional partons starting from the scale p\^- There 
are now two cases 


(a) If n = the trial emission is always kept and the dipole cascade is allowed to 
continue down to the cutoff p^^ and the event is accepted. 

(b) If n < N, and all parton pairs pass the cut, Qq, used in the MEG, the event is 
rejected and we restart from step 1. If any of the partons fail the cut, the trial 
emission is accepted and the dipole cascade is allowed to continue down to the 
cutoff p^^ and the event is accepted. 


5. Results 


To test our algorithm, we have generated W'’'-|-njet events, with n < N = 4 with the 
MadGraph/MADEvent program [49] for a pp collider at a total energy of 1960 GeV, 
ie. corresponding to the Tevatron run II. The longitudinally invariant fc_L-algorithm was 
used^ to regularize the cross section, using cutoffs Ecut = 12,17 and 22 GeV. We used 
the CTEQ6L [50] PDF parameterization using as scale. E^^^ was also used as the 
scale in cts- The event was generated with unit weight and was then reweighted according 


to the algorithm in section 4.3. To avoid wildly fluctuating weights due to the ratios of 


PDFs, arising from situations in which the events in MadGraph had large x values where 
the PDF is very small, the W"*" was required to have limited rapidity, |p| < 2.5, and all 
partons were required to have a limited pseudorapidity, [p] < 2.5. For the same reason, 
when constructing the remnants, both the valence- and sea-quark alternatives were used 


^Using MDDE=4211 in the KTClus program [47]. 
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(a) 


(b) 





Figure 3: Cross section as a function of transverse momentum of the W for standard AriaDNE 
compared to the first order matrix element correction for the different cutoffs (a) 12 GeV, (b) 
17 GeV and (c) 22 GeV. Plot (d) shows the ratio between the distributions including matrix 
element corrections and standard AriaDNE. 


with appropriate weights, rather than choosing between the two as described in section 

3l. 

In the following we will use the notation MEA^PS for results from our new algorithm 
using W'’'+0,... ,W'’“+Af jets from MadGraph (the individual contributions from + n 
jets are denoted ME’^A^PS), while Ariadne will denote results obtained with the default 
Ariadne treatment. 

As mention above, Ariadne by default already has a matrix element correction for 
the first emission in W production^. Hence, as a first test of our new algorithm is to run it 
with = 1, in which case we should get the same result as the standard default Ariadne. 
In fact, had the construction procedure been exact, the results would be exactly the same. 
Of course, the construction can never be really exact, but for only one additional jet it is 
fairly close. 

In figure ^ we show the p± spectrum of the for the new algorithm with A^ = 1, 

®In fact we discovered a small bug in the default Ariadne treatment, related to a mismatch between 
the invariant pxand the actual transverse momentum in gluon emissions. This bug has been fixed in the 
latest release of Ariadne. 
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(a) 


(b) 




Figure 4: Cross section as a function of transverse W momentum for standard AriaDNE, AriaDNE 
with fourth order matrix element corrections and Pythia with 1 GeV^ and 4 GeV^ intrinsic k± 


MEIPS, compared to default Ariadne. Clearly the agreement is very good. In particular 
we note that there is no significant discontinuity or other strange behaviors around the 
different cutoffs used. This agreement is not trivial, since the MEIPS curve is a sum of 
W^+0 and W^+1 jet event from MadGraph. The fact that there is a small contribution 
of W^+0 event above the cutoff (and vice versa, some W'’'+l events below the cutoff) is an 
effect of the added cascade. The effect of the cutoffs is not completely invisible, however, 
as is clear in figure id, where we enhance the effect by showing the ratio of p^-spectra of 
MEIPS and Ariadne for the different cutoffs. 

We can now proceed with some confidence to investigate our new algorithm also for 
higher parton multiplicities in the MEG. Here we will, of course, expect differences w.r.t. 
Ariadne. These differences would be the ones desired from replacement of the products of 
splitting functions with the exact tree-level matrix elements, and also from the additional 
processes not present in the dipole cascade, such as the one in figure For small scales, 
these differences should be small and we would still like to have a smooth behavior of 
any observable sensitive to the cutoff used in MadGraph, at least for small enough Ecut- 
There may also be differences arising from deficiencies in our algorithm, since there are 
now additional construction steps possible. 

In figure ^a we again show the W"*" p±_ spectrum, but now using ME4PS and comparing 
with Ariadne and also with the default PYTHIA parton shower (which also includes a tree- 
level matrix element correction for the first emission). We find that there is now an increase 
at large p±_ for the ME4PS case, which is attributed to the desired higher-order effects. We 
note that there is still no dramatic discontinuity around the cutoff. 

There is distinct difference in the small-px behavior for the PYTHIA distribution. The 
peak in PYTHIA is shifted towards smaller as compared to Ariadne, with or without 
the new matching algorithm. This is a known problem with PYTHIA, which need an 
uncomfortably large intrinsic transverse momentum of the proton to reproduce data. This 
is shown in more detail in figure ^b where we focus on the small-p^ part of the spectrum 
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(a) 


(b) 






Figure 5: Cross section as a function of the scale where the jets merge in the fcj_-algorithm, where 
di denotes the scale when i jets merge to * — 1 jets. Standard AriaDNE is compared to fourth order 
matrix element corrected distributions at a 12 GeV cutoff 

and where we have two curves for Pythia, one with an average squared intrinsic transverse 
momentum, of 1 GeV^ and one with 4 GeV^. Note that also Ariadne has an intrinsic 

transverse momentum, but this is at a typical non-perturbative value of = 0.36 GeV^, 
but the increased possibility of radiating gluons in the dipole cascade, especially in the 
direction of the remnants, will give the slightly harder spectrum. From reference [51] we 
know that PYTHIA can only describe data with the higher intrinsic transverse momentum^, 
which is quite close to Ariadne ® There is still a difference in shape and with the increased 
statistics collected in Tevatron Run II, it may be possible to distinguish between the two. 

Next we want to check the cutoff sensitivity also for higher jet rates. We do that 
by taking the final events on parton level and cluster them with the same /c^-algorithm 
which was used for the regularization in MadGraph, and then look at what value of the 
resolution variable, dn, an event is clustered from n-jets to n — 1-jets. In figure we show 
such distributions for n = 1,2,3 and 4 for Scut = 12 GeV. We also show the individual 
contributions from different parton multiplicities delivered by MadGraph. We see that 

^In later Pythia releases the higher value is the default. 

®We have not compared directly with data here due to uncertainties about the corrections made to the 
data. 
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(a) 


(b) 




(c) (d) 




Figure 6: ME corrected distributions and Pythia divided by AriaDNE as a function of the jet 
merging scale in the fcj^-algorithm for all different cutoffs used. 


there is a clear difference between ME4PS and Ariadne for large dn values, which is 
expected from the improved treatment of events with several hard jets. We note that there 
is a rather smooth transition across Ecut- In figure we also show results for Ecut = 17 
and 22 GeV, now presented as ratios between ME4PS and Ariadne. As expected we here 
see more clearly where the new matrix element treatment sets in above Ecut • In figure ^ we 
also show the ratio between PYTHIA and Ariadne and we find that, as compared to the 
effects of the matrix element corrections, the difference between the two cascades is small. 

Next we want to see if the features of the Ariadne resummation are reflected in our 
new algorithm. As noted before, the emission of gluons are allowed in a larger phase 
space region in Ariadne as compared to a conventional PSEG, hence the no-emission 
probabilities should be affected. Also, in a conventional DGLAP-based initial-state PSEG, 
the parton closest to the W is also the hardest one. This is not the case for Ariadne where, 
effectively, contributions of emissions with lower p± between the hardest parton and the W 
is taken into account, as illustrated in figure ^a. One observable which may be sensitive to 
this difference is the pseudorapidity difference between the W and a jet, Arywj, which then 
should be enhanced for large Ar/wj in Ariadne as compared to PYTHIA. This is also the 
case as shown in figure ^b. Of course, the exact tree-level matrix element will also contain 
contributions such as the one in |^, and we see that the enhancement at large Ar/wj is even 
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Figure 7: Part (a) shows a typical diagram contributing to the pseudorapidity difference between 
the W and the hardest jet. Diagrams with soft jets between the W and the hardest jet can not 
be generated through DGLAP evolution, but are included in the matrix element corrections. Plot 
(b) is a normalized distribution of the difference in pseudorapidity between the hardest jet and the 
W, where the jet is defined using the fc_L-algorithm with a 12 GeV cutoff and the hardest jet has 
a transverse momentum greater than 40 GeV. The distribution is shown for AriaDNE (full line), 
Pythia (dotted line), and ME4PS with a 12 GeV and 22 GeV cutoff (long- and short-dashed lines 
respectively). 


more significant, when ME corrections are added to Ariadne. It would be interesting to 
see if including CKKW corrections also in PYTHIA would bring it closer to Ariadne and 
ME4PS. This could be expected since including higher order corrections could in the end 
make things more insensitive to the particular kind of resummation used. 

One of the advantages of the matrix element corrections is that correlations between 
hard partons are more accurately described. This may be important when eg. estimating 
backgrounds to different searches. We will here consider the background to top production 
at the Tevatron for the semi-leptonic channel which corresponds to W-|-4-jets. In a realistic 
top search one would use identified b-jets, but since our MadGraph events do not include 
b-quarks we look at W-|-4-jets in general. 

In figure ^ we show the W-|-4-jets background to the top-mass distribution. We ob¬ 
tained it by using the /c^-algorithm to cluster four jets and required that the jet scale 
was above 12 GeV. Form these we found the two jets ji and with an invariant mass 
mi 2 closest to the W mass. If no jets with |mi 2 — mw| < 20 GeV were found, the event 
was rejected. Then we selected a third jet, j 3 , so that the difference |mi 23 — mw 4 | was 
minimized, where mi 23 is the invariant mass of jets 1, 2 and 3, and mw 4 is the invari¬ 
ant mass of the W^and jet 4. If the difference |mi 23 — mw 4 | < 20 GeV then the event 
was accepted and the constructed top mass is defined as the average of mi 23 and mw 4 - 
Clearly the total cross section would be underestimated by the leading-order predictions of 
Pythia and Ariadne. In figure we therefore only show the normalized shape and find 
that Pythia and Ariadne are quite similar and that no significant change is introduced 
by the matrix element correction to Ariadne. We also show the result from using the 
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Figure 8: This is the background for a search for the top quark. A normalized distribution as 
a function of the mass of the top quark for AriaDNE (full line), Pythia (short-dashed line), and 
ME4PS with 12 GeV cutoff (long-dashed with error bars). Also shown is the results from the pure 
tree-level matrix element without parton showers added using a 12 GeV cutoff (dotted line). 



^12 

Figure 9: Normalized distribution of the azimuthal angle difference between the two hardest jets 
defined using the A:_L-algorithm with a 12 GeV cutoff. The distribution is shown with AriaDNE (full 
line), Pythia (short-dashed line), and ME4PS with 12 GeV cutoff (dashed line with error bars). 


tree-level 4-jet matrix elements directly, without reweighting and adding a cascade, and 
find no large difference with the parton shower approaches. Possibly the fall-off with mt 
can be said to be somewhat weaker for the pure matrix elements. 

To focus more specifically on angular correlations, we finally look at the azimuthal angle 
between the two hardest jets, (/)i 2 . This observable is important for understanding how 
higher order emissions influence the transverse momentum of the W, 4’i2 ~ t. 
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the emission of a second jet decreases the p_lW) while for (j)i 2 ~ 0 the p±w is increased. 
In Ariadne, after a first emission of a gluon, a second gluon will be radiated isotropically 
in azimuth in the rest system of the radiating dipole. Since this dipole is boosted in the 
direction of the first gluon, we expect that the second gluon is more likely to go in the 
same direction. This is not true for Pythia, where successive initial-state emissions are 
uncorrelated in azimuth (with the recoils, the net effect is a bias towards large <^ 12 )- In 
figure ^ we see the (normalized) (j)i 2 distributions, and indeed we find that Ariadne is 
more biased towards ^12 ~ 0. Adding matrix element corrections removes this bias and 
brings the distribution closer to the PYTHIA result. Hence this indicates that the azimuthal 
correlations in standard Ariadne are not very well modeled. 

6. Conclusions 

We have presented a way to implement a CKKW procedure for combining events gener¬ 
ated according to tree-level matrix elements for W-production with the dipole cascade of 
Ariadne. Although the basic principles are fairly simple, the details of our procedure is 
rather involved, which is mandated by our aim to become as insensitive as possible to the 
cutoff needed in the matrix element generation. 

Our strategy is to take any partonic state generated by a MEG and try to find a 
likely history of emissions which Ariadne would have performed in order to generate this 
state. Rather than just using the constructed emission scales to calculate analytic Sudakov 
form factors to reweight the states, as is done in the original CKKW procedure, we find 
exactly the Sudakov form factors Ariadne would have used. In addition we reweight the 
states with the parton densities functions and the so-called soft suppression function which 
Ariadne would have used. 

The PDF reweighting means that the overall normalization of cross sections are still 
given by the leading order diagrams used by standard Ariadne. However, we expect 
much improvement of the shapes of final state distributions as compared to the standard 
parton shower description. In [7] it was found that the CKKW procedure reproduces well 
the shapes of distributions obtained by a NLO program such as MC@NLO [52], although 
the overall normalization needed to be adjusted with a K-factor. Although not explicitly 
checked in this report, we expect that this will also hold for our implementation. 

We have presented several investigations into how the Ariadne program is improved by 
adding matrix element corrections. In some cases we also compared to the PYTHIA parton 
shower to get some insight into how well these standard cascade programs reproduces 
higher order matrix elements. 

In one case we looked at the azimuthal correlation between the two hardest jets, and 
found that the difference between Ariadne and PYTHIA was large. When corrected with 
matrix elements, Ariadne came much closer to PYTHIA, indicating that such azimuthal 
correlations are not handled very well in standard Ariadne. 

In our quasi-realistic top-background observable we found that Ariadne and PYTHIA 
were quite close and that no drastic effect was obtained by including matrix element cor¬ 
rections. 
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For the W-jet rapidity correlation we again found clear differences between Ariadne 
and Pythia, and that these were even enhanced when correcting Ariadne with matrix 
elements. This indicates that such correlations are not very well described by PYTHIA, 
while standard Ariadne does a better job, although it can be improved. 

We believe that the rapidity correlations indicate that non-ordered evolution is of 
importance for W-production at the Tevatron. Such evolution is expected to be important 
in small-x processes, and the fact that it shows up here, where x ~ mwl\fS K, 0.04, may 
be somewhat surprising. We also believe that the inclusion of non-ordered evolution is 
why Ariadne is able to reproduce experimental data on the small-p_L-distribution of the 
W and TP, distributions which can only be described by PYTHIA if an uncomfortably large 
intrinsic transverse momentum is added. 

The fact that matrix element corrections can give us hints about where unordered 
evolutions may become important, is an indication that it would be very interesting to 
implement CKKW also for DIS® and compare with HERA data. Also, at the LHC where 
W-production may be argued to be a true small-x process (x mwiVS P, 0.006), it should 
be interesting to study matrix element corrections. In fact also Higgs production at the 
LHC may be considered to be a small-x process. We will come back to these processes in 
future publications. 
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Figure 10: Symbolic pictures of possible steps which are possible when constructing possible inter¬ 
mediate partonic states from events generated by a MEG. The different possibilities are described 
in the text. The dashed lines indicate colour-connections. 
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A. Appendix: Construction 

Here we describe in some detail the different kinds of steps possible when constructing 


intermediate partonic states from events generated by a MEG. In figure 1^ the different 
steps are shown schematically. 

(a) Pi—g 2 —P 3 —> Pi—P 3 - A gluon, colour-connected to two non-remnant partons, pi and 
P 3 , is constructed to a single dipole between pi and p^. The splitting function is given 
by one of eqs. ( ^ ) - ( |2.11 ) depending on whether pi and ps are gluons or quarks. 
The scale is given by the invariant p± in eq. (^^). No PDF ratio is relevant. In the 
rest frame of the construction, pi will retain its direction if it is a gluon and pa is a 
quark, and vice versa. If both pi and ps are gluons or both are quarks, the direction 
of Pi is rotated away from the original pa direction with an angle /3x^/(x^ + 3:|), 
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where /3 is the original angle between pi and pa. This corresponds to the standard 
recoil treatment for gluon emission in AriADNE. 


(b) pi—q 2 qs —> Pi— g- A qq-pair is constructed into a gluon as long as they are each 
others anti-particles, are connected to different strings, and the end-points of these 
strings are not remnants of the same incoming hadron. The splitting function is given 
by eq. (2.13). The scale is given by the invariant p_L in eq. ( |2.6| ). No PDF ratio is 
relevant. In the rest system of the construction, pi will retain its direction. 


Pi-g2-rz 


Pi—r^: A gluon, colour-connected to one remnant, and one non¬ 


remnant parton, pi, is constructed to a single dipole between pi and rs. The splitting 
function is given by eqs. ( |2.g| ) or ( 2.10| ) depending on whether pi is a quark or a gluon. 
If r 3 is one of two remnants of the same hadron (this corresponds to an extracted 
gluon), the PDF ratio is taken to be 0, otherwise it is 0/z. If a W is present in 
the event and it is close to g 2 , the transverse momentum of the gluon in the center 
of mass system of pi and is given to the W, and the longitudinal momentum is 
absorbed by pi and r^. The scale is given by the invariant p± in eq. (|2.6|) (calculated 
as if no W was close, ie. the transverse momenta of the gluon is transfered to pi and 
rs with the weight x\l{x\ -|-x|)). Here close means that p+g < p+w and p_g < p-W) 
where p±w is calculated for the constructed W momenta. If there is no W close by 
the momentum of the gluon is shared by pi and r^, where retains its direction. 
The scale is given by the invariant p_L. 


(d) ri—g 2 —ra —> ri—r^: A gluon connected to two remnants, one from each incoming 
hadron, is constructed to a single dipole between the remnants. The splitting function 
is given by eqs. (2.£). When calculating the scale a fraction xfj{x1 + x|) of the 
transverse momenta from the gluon is transfered to each of the remnants and the 
scale is given by the invariant pj_ in eq. ( |2.6D . The PDF ratio is given by the product 
of the 0 on each side, divided by z if the corresponding remnant is not one of two 
remnants of the same hadron. The transverse momentum of the gluon is transfered to 
the W if one is present, otherwise it is transfered to the hard subsystem containing 
the rest of the non-remnant partons in the event. The longitudinal momentum is 
divided between ri and r^. 


(e) Hi Tg q 2—^3 —> Hi r'^ h^: This corresponds to the inverse of an initial-sate g ^ q 
splitting. For a quark, q 2 , connected to a remnant, r^, and a hard subsystem. Hi 
(which contains the W if present), connected to another remnant, rg, from the same 
incoming hadron and arising from the extraction of a corresponding anti-quark, q^, 
a hadron, / 13 , is formed from q 2 and r^. The splitting function is the standard 
AltarellDParisi one, Pg^q(z). The scale is the squared transverse momentum of q 2 
in the rest frame of the event. The PDF ratio is the same as would have been used 
in a conventional parton shower. The transverse momentum of q 2 is transfered to 
Hi, and the longitudinal momentum is shared between Hi, rg and / 13 . The relative 
sharing of longitudinal momenta between rg and is the same as for the original rg 
and r 3 . 


(f) Hi q 2 r 3 (/ 13 ) —> Hi r^: A quark, q 2 , which may have been extracted from a 
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hadron resulting in a remnant may be absorbed into a the remnant, constructing 
an initial-state q ^ g splitting. The remnant is split into two, possibly together with 
a remnant hadron, h^, if q2 was a sea-quark. The splitting function is the standard 
Altarelli-Parisi one, Pq^g(2;). The scale is the squared transverse momentum of q2 
in the rest frame of the event. The PDF ratio is the same as would have been used in 
a conventional parton shower. The transverse momentum of q2 is transfered to the 
spectator hard subsystem, Hi, and the longitudinal momentum is shared between 
Hi, rs and r'^. The relative sharing of longitudinal momenta between ra and r'^ is the 
same as for the original rs and /13 if was present, otherwise the momenta is shared 
as is normally done in Ariadne when a gluon is extracted from a hadron. Note that 
there is no corresponding emission in AriADNE. 



